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Abstract 

We study the intrinsic properties of attractors in the Boolean dynamics in complex 
network with scale-free topology, comparing with those of the so-called random 
Kauffman networks. We have numerically investigated the frozen and relevant nodes 
for each attractor, and the robustness of the attractors to the perturbation that flips 
the state of a single node of attractors in the relatively small network (N = 30 ~ 
200). It is shown that the rate of frozen nodes in the complex networks with scale- 
free topology is larger than that in the random Kauffman model. Furthermore, we 
have found that in the complex scale-free networks with fluctuations of in-degree 
number the attractors are more sensitive to the state flip of a highly connected node 
than to the state flip of a less connected node. 
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1 Introduction 



Dynamics of Boolean networks is often used as a model for genetic networks 
inside cells, in which the genetic states are represented in terms of the language 
of attractors [1,2]. A Boolean network consists of iV nodes, each of which 
receives k{ inputs such that the degree is ki. In the so-called Kauffman model 
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- a random Boolean network (RBN) model, each node receives a certain fixed 
number of inputs such that the degree is ki = K. 

However, in real systems each node has a different number of inputs. As a 
more realistic modeling of biological systems we expect that the fluctuation 
of the number of input-degree is treated as a random variable with a proba- 
bility distribution function such as the inverse power-law distribution or the 
exponential distribution or the Poisson distribution. In fact, the in-degree dis- 
tribution appears to be exponential in E.coil and to be power-law in yeast 
[3-9]. 

In the present report, the number of in-degree fcj at the %— th node is deter- 
mined by the preferential attachment rule [10] that makes the system a com- 
plex network with scale- free topology (SFRBN). The details of the method for 
generating the networks have been given in our previous paper [11]. It should 
be stressed that in the present paper we do not deal with large networks that 
the difference between the function forms of the distributions becomes clear 
numerically. The network size that we deal with is at most N = 30 ~ 200, 
focusing on the intrinsic properties of Boolean dynamics in the complex net- 
works. 



2 Model 



The initial values for the nodes are chosen randomly and are synchronously 
updated in the time steps, according to the connectivity {hi} and the Boolean 
functions {/j} assigned for each node in the network as, 

(T i (t + l) = Ma il (t),a ia (t),---,(7 iki (t)), (1) 

where % = 1, N and <7; G {0, 1} is the binary state. All trajectories starting 
at any initial state run into a certain number of attractors(i.e. points or cycles). 
We study the directed RBNs, the directed SFRBNs, and the directed SFRBNs 
throughout this paper. 

Figure 1 shows the time-dependence of the state which constitutes an typical 
attractor in the SFRBN with (k) = 2. 

Fig. 2(a) shows the histogram of the length £ c of the attractors in the RBN 
with K = 2 and that in the SFRBN where the average degree of nodes (k) = 2. 
Fig. 2(b) shows the median value fa of the distribution of state cycle lengths 
with respect to N, the total number of nodes. Apparently, the distribution 
of the attractor lengths in the SFRBN is much wider than that in the RBN, 
and the attractor length has longer period than that in the RBN. This is 



2 



50 



CD 
-Q 

E 

c 

CD 
T3 
O 

c 









■III 1 . ■ 1 ¥ J 











time step 



88 



Fig. 1. (Color online) Space-time diagram for an typical attractor with £ c = 88. in 
the SFRBN with {k) = 2. The vertical axis denotes the node number that is in the 
high connectivity order. The dense or gray color corresponds to the binary values 
of the state. The static nodes of the attractor are frozen nodes. The network size is 
N = 50. 
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Fig. 2. (Color online) (a) Histogram of the length £ c of state cycles is shown for 
the RBNs and the SFRBNs, respectively, where the network size is N = 200. Each 
histogram is generated by 10 3 different sets of the Boolean functions and five differ- 
ent network structures. The maximum iteration number of the Boolean dynamics is 
10 5 until the convergence to the cycle is realized, (b) Semi-log plots of the median 
value fh of 1000 samples of the lengths of the state cycles with respect to the total 
number N of nodes for the directed RBNs and SFRBNs. 
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directly related to the diversity of attractors in the SFRBNs, which is of great 
importance for the stability of living cells. We investigated the function form 
fh(N) in more detail in the previous paper [11], and found that the function 
form fh(N) asymptotically changes from the algebraic type m(N) oc N a to 
the exponential one as the average degree (k) goes to (k) =2. 

Here we have a question: How does the characteristics of the attractors depend 
on the network topology? In the present paper we study the difference between 
attractors in the RBN and the SFRBN without a bias, focusing on the frozen 
nodes and the robustness of attractors against the external perturbations. 



3 Frozen Nodes of Attractors 

In this section, we investigate the so-called frozen nodes of attractors whose 
values remain constant through a given trajectory of the attractors [1]. Frozen 
nodes arise through canalizing the Boolean functions and the homogeneity 
bias. 

We count the number of frozen nodes Nf for each attractor and plot the 
histograms for some cases in Fig.3. The remarkably different peak structure 
exists between the cases in the SFRBN of (k) = 2 and the RBN of K = 2. 
The distributions in the SFRBN have a peak around Nf ~ N/2, while the 
distributions in the RBN are broad with a peak at Nf = N. Note that Nf = N 
corresponds to the point attractors that all nodes are frozen. 

The fact that Nf is relatively smaller in the SFRBN would make the attractor 
period larger than that in the RBN, as seen in Fig.2. However, role of the 
frozen nodes is not so clear in the network because all frozen nodes also are 
connected to the network and might influence on the attractors. In the next 
section, we investigate the significance of each node in the respective attractor. 



4 Robustness of Attractors to State Flip 

One of the important properties in the scale-free topology is the existence of 
the highly connected hub node as seen in the yeast synthetic network and so 
on. In this section, we investigate the robustness of attractors to an external 
perturbation caused by an inversion of the binary state of a single node. We 
consider an attractor of period £ c and flip the state of the single node at time 
t G [1, £ c ] as a perturbation. The perturbation to the trajectory of the attractor 
may leap from the trajectory of the original attractor to another one, i.e. the 
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Fig. 3. (a) Histograms of the number of frozen nodes Nf for 1000 attractors in RBNs 
with K = 2 and SFRBNs with (k) = 2. The network size is (a) TV = 50, (b)iV = 100 
and (c)N = 200. The scale out data at Nf = N are not shown in the figures. 



attractor shift. The high homeostatic stability implies low reachability among 
different attractors. 



We investigate the probability R s that the attractor remains in the original 
attractor under the inversion of the single node state [2,14], which is the rate 
returning to the original attractor under the inversion. Here we call the rate 
the robustness of the attractor. 



Figure 4 shows the robustness R s of attractors with £ c = 55 in the SFRBN and 
£ c = 48 in the RBN, respectively. It follows that in the SFRBN the number 
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Fig. 4. The rate R s of returning to the original attractor as a function of numbered 
nodes in the order of the number of in-degree k{. We used N = 100 in both the 
SFRBN with (k) = 2 and the RBN with K = 2. The periods of the attractors are 
(a)4 = 9 in the SFRBN and £ c = 9 in the RBN, (b)4 = 55 in the SFRBN and 
£ c = 48 in the RBN, respectively. The horizontal axis "Node number" denotes the 
node number in the order of the number of input-degree. 

of "active nodes" (R s < 1) is much more larger than that in the RBN. On 
the other hand, in the RBN the perturbation to the active nodes influences 
effectively the shift of attractor (R s < 0.6) although the number of active 
nodes is not so many [13]. As a result, in the SFRBN the perturbation to the 
highly connected hubs may give rise to the attractor shift, comparing with the 
one to the less connected nodes. 

The robustness of some attractors in the SFRBN with (k) = 4 and in the 
RBN with K = 4 is given in Fig. 5. Although the whole structure is almost 
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Fig. 5. Robustness R s for some attractors of the SFRBN (k) = 4 and the RBN with 
K = 4 in the network size of N = 30. The periods of attractors are l c = 73 in the 
SFRBN and £ e = 62 in the RBN. 

similar to that for the cases in Fig. 4, it is found that the effect of inversion of 
the single site state on the attractor shift is relatively small compared to the 
cases of (k) = 2 and K = 2. There is a tendency that the attractors become 
more robust to the perturbation as the average number of input-degree (k) 
increases. 



5 Summary and Discussion 



In summary, we have studied the Boolean dynamics of the Kauffman model 
with the directed SFRBN, comparing with the ones with the directed RBN for 
the relatively small network size. In this study we investigated some intrinsic 
properties of attractors between the RBNs and the SFRBNs, focusing on the 
frozen nodes and the robustness to a perturbation. The obtained results are 
as follows. 

(i) The number of frozen nodes in the SFRBN is smaller than that in the 
RBN and the property reflects on the much more widely distributed attractor 
lengths. 

(ii) The perturbation to the highly connected hubs may give rise to the at- 
tractor shift in comparison to the less connected nodes. 

(hi) The attractors becomes more robust to the perturbation as the average 
number of input degree (k) increases. 

Although in this report we did not show the details of the attractor shifts by 
the perturbation, we will present the details of the numerical results for the 
diagram of transition among the attractors and the robustness to perturbation 
in our forthcoming paper [13]. 
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Robustness against genetic mutations and environmental perturbations is one 
of the universal features of biological systems. And the robustness is important 
for understanding evolutionary processes and homeostasis of gene regulatory 
networks [17-19]. However, in this report we investigated only robustness to 
the single site inversion. For the purpose of the study, the other robustness of 
the SFRBN might become significant; for instance, robustness of attractors to 
the change of Boolean functions and to the breakdown of the network struc- 
tures. We expect that such a study on the robustness of attractors provide 
some insights into important biological phenomena such as cellular homeosta- 
sis and apotosis. 

Actually Aldana et al investigated the small SFRBN (N ~ 15 — 20) and found 
that the robustness of the ordered phase to the network damage is lower than 
that in the RBN [14] The result implies that there exists a possibility of evo- 
lution though mutation even in the ordered phase, despite of the Kauffman's 
conjecture that life evolves in "edge of chaos". 

Moreover, recently, an interesting network model, the so called feedback net- 
work has been proposed by White et al [20]. The feedback networks can be a 
good model for describing the autocatalytic chemical reactions and the kin- 
ship, and so on, because the node selection, the search distance and the search 
path of networks are controlled by the attachment, the distance decay and the 
cycle formation parameters. It is interesting to investigate the features of the 
Boolean dynamics in the feedback networks from the point of view of frozen 
nodes and robustness [13]. 
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